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^ ' Abstract 

The integrable structure, recently revealed in some classical problems of the 
theory of functions in one complex variable, is discussed. Given a simply connected 
domain in the complex plane, bounded by a simple analytic curve, we consider the 
conformal mapping problem, the Dirichlet boundary problem, and to the 2D inverse 
potential problem associated with the domain. A remarkable family of real-valued 
functionals on the space of such domains is constructed. Regarded as a function 
Q> I of infinitely many variables, which are properly defined moments of the domain, 

' any functional from the family gives a formal solution to the problems listed above. 

■ These functions are shown to obey an infinite set of dispersionless Hirota equations, 

i This means that they are r-functions of an integrable hierarchy. The hierarchy is 

identified with the dispersionless limit of the 2D Toda chain. In addition to our 
previous studies, we show that, with a more general definition of the moments, this 
connection is not specific for any particular solution to the Hirota equations but 
reflects the structure of the hierarchy itself. 
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>^ ■ 1 Introduction 



Recently, it was found [Q, ^ ^ that some problems of complex analysis, such as conformal 
mapping problem, Dirichlet boundary problem, and 2D inverse potential problem, have 
a hidden integrable structure. It turns out that there is an infinite hierarchy of non-linear 
integrable equations such that variations of sought-for solutions to the problems hsted 
above under deformation of given data can be found by solving the hierarchy. The latter 
is a certain scaling limit, usually referred to as dispersionless limit, of the 2D Toda chain 
hierarchy. 

The conformal mapping problem is to find a conformal map from a given simply 
connected domain in the complex plane to the unit disk. By the Riemann mapping 
theorem, such a map does exist and, after imposing some simple normalization conditions, 
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is unique (see, e.g., 0). To be specific, let 7 be a closed analytic curve in the complex 
plane and D_ be respectively the interior and exterior domains with respect to the 
curve. By w{z) we denote the conformal map of D_ to the exterior of the unit circle 
normalized in such a way that 00 is taken to 00 and the derivative at 00 is real positive. 

The Dirichlet boundary problem in sounds as follows: given a function il){z) on 
the curve 7, to find a harmonic function f{z, z) in D_ such that /(z, z) = ip{z) for z G 7. 
The solution is given by the formula 

f{z,z) = --i^{Od^G{zX)dC (1) 

where G{z, () is the Green function of the Dirichlet problem in D_. The Green function 
is uniquely characterized by the properties 0: (a) It is regular at infinity and harmonic 
everywhere in D_ but at the point z = ( where it has the logarithmic singularity: 
G{z, = log 1^; — CI + 0(1) as 2; — > C) (b) It is equal to zero on the boundary dD_ = 7, 
i.e., G{z, = for all ( if z ej. 

To specify what we mean by the inverse potential problem, we suppose, in the same 
setting as above, that the domain is filled with a continuously distributed background 
electric charge with uniform density. As is proved in the potential theory, the potential, 
$(2;, z), generated by the charge, and its first derivatives (components of the electric field) 
are continuous functions at the boundary, i.e., = $~, 9^$'*' = dz^~ , 267. Here 
$^ is the function $ restricted to D±. $^ is harmonic everywhere in but at infinity, 
where it has a logarithmic singularity, while is harmonic up to a term proportional 
to l^rp. Therefore, one can represent by Taylor series (multipole expansions) in a 
neighbourhood of and 00 respectively. We suppose, without loss of generality, that 
is in The coefficients of the multipole expansions are harmonic moments of D_ and 
D^, respectively. The inverse potential problem is to restore the shape of the domain 
given one of the series, either or $~. For definitness, we assume that the coefficients 
of are given and address the problem how to find the coefficients of $~ . (If one knows 
both sets of coefficients, the curve can be easily reconstructed.) For a review of inverse 
problems of potential theory see, e.g., and references therein. 

The three problems are closely related to each other. Given, say, the conformal map 
wlz), the Green function can be found explicitly via the formula 



G{z,C) = log 



w{z) — w{() 



w{z)w{() - 1 



(2) 



where the bar means complex conjugation. In turn, the holomorphic in z part of the 
Green function is the conformal map, though normalized in a different way. To obtain 
the map w{z), one should take the holomorphic part of G{z, 00). Concerning the inverse 
potential problem, we note that it is the Green function that provides a relation between 
and . Indeed, let us modify the potential function and introduce $ which is $ with 
subtracted logarithmic singularity at infinity, so that $^ is harmonic in D_ including 00. 
Then, according to (|^), we can write 

z) = --i $+(C, C)dcG{z, QdC (3) 
m J J 
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Let us summarize the main results of |jl], ^, ^. Let t^, k > l,he the harmonic moments 
of D_, and to be the area of (divided by vr): = — ^ z^'^dPz, = ^ /d+ d'^^- 
Note that are in general complex numbers while to is real. In the subsequent formulas, 
the complex conjugate moments, ik-, are considered as independent variables. Introduce 
the operators 

D{z) = E ^5*. , m = y: (4) 
k>i k>i 

and 

V{z,z) = dt, + Diz)+D{z) (5) 

There exists a real-valued function F of to, {tk}, {tk} such that the conformal map w{z) 
is given by 

w{z) = zexp [—dlF - dt,D{z)F^ (6) 

This function admits an explicit integral representation. Moreover, the Green function 
of the Dirichlet boundary problem is represented as 

G{zi,Z2) =\og\z^^ - Z2^\ + h:){z^,z^)V{z2,Z2)F (7) 

Let us note that a formula of this type was first conjectured by L.Takhtajan. See for 
a rigorous proof and discussion. 

The function F is the dispersionless limit of the logarithm of the r-function of the 
2D Toda chain hierarchy. It obeys the dispersionless Hirota equations 

ZiZ2 (l - e--°(^i)^(^2)F^ ^ ^at^{dto+D{z^)+D{-Z2))F ^^g^| 

The dispersionless Toda hierarchy is an example of the universal Whitham hierarchy 
introduced in It is a multi-dimensional extension of the hierarchies of hydrodynamic 
type . Equations of the Toda hierarchy are known to have infinitely many solutions. 



The solutions can be parametrized |Tl[ by canonical transformations in a two-dimensional 
phase space in such a way that any solution corresponds to a canonical pair of functions. 
(In the canonical pair is called twistor data of the solution.) In fact, this is equivalent 



to the characterization of the solutions via string equations. 

In the approach developed in 0, it was one particular solution to the dispersionless 
2D Toda hierarchy that described conformal maps. The place, if any, of other solutions in 
the conformal map context was not clear. In this paper we address exactly this question. 
Extending the previous result, we show that the same class of conformal maps can be 
described by any generic solution to the hierarchy. 

This description is achieved by another, more general definition of moments of the 
domain. The moments are again the independent flows of the hierarchy. However, their 
relation to the shape of the domain is different, and some additional data is required for 
their definition. The additional data is a real analytic function on the complex plane, 
a{z,z). In the context of the inverse potential problem, this function is density of the 
background charge, which was constant in the standard formulation of the problem. More 
precisely, we suppose that a{z, z) is given in the whole complex plane and fixed once and 
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for all while the shape of the domain may vary. The background charge has density 
a{z, z) inside the domain and outside. The inverse potential problem is posed in the 
same way as before. The moments, however, are now defined with respect to the density 
cr. In the sequel we refer to this problem as generalized inverse potential problem. 

The result shows that it is the hierarchy of Hirota equations (rather than particu- 
lar string equations or other solution-specific ingredients) which reflects the integrable 
structure of conformal maps and the Dirichlet problem. 

In Sec. 2 we introduce the function F with the help of a variational principle and find 
its first derivatives. The first derivatives give a formal solution to the generalized inverse 
potential problem. Formal solutions to the conformal map problem and the Dirichlet 
problem are given by second derivatives of F. There is an infinite set of relations for 
the second derivatives, which can be combined into the (dispersionless limit of) Hirota 
equation. This stuff is discussed in Sec. 3. At last. Sec. 4 is devoted to an alternative 
approach to the dispersionless integrable hierarchy, which allows one to prove the Lax- 
Sato equations for conformal maps. 

2 Generalized inverse potential problem 

The inverse potential problem is perhaps the most convenient context to introduce the 
r-function for curves. Following the idea of [Q, we define the r-function with the help 
of an electrostatic variational principle. The set-up and the notation are the same as in 
Sec. 1. 

We recall that the domain Dj^ is filled with a background electric charge. Let us 
consider non-uniform charge distributions such that density of the charge is characterized 
by a real analytic function a{z,z). We assume that this function is defined all over the 
complex plane, the charge density being equal to —a inside D+ and outside. With a 
fixed function o", the 2D electrostatic potential $ generated by the domain is a functional 
of the shape of the domain. The potential $ obeys the Poisson equation 




if z = X + iy G D. 
if z = X + iy E D 



+ 



(10) 



The potential $ can be written as an integral over the domain D 




(11) 



In D_, the potential is a harmonic function represented by the series 



$-(^, z) = -2tolog|z| + 27^ey 



(12) 




TT JD+ 

are moments of the interior domain Dj^ with respect to the density cr, and 



(13) 




(14) 
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is the total charge of D^. 

In the potential (p!TD is represented by a series which has both harmonic 
and non-harmonic terms. The latter is due to the background charge. We denote the 
non-harmonic part by U{z, z). The function U obeys the Poisson equation dzdzU{z, z) = 
a{z, z) and is fixed by the condition that its expansion around (recall that G -D+) 
does not contain harmonic terms. For example, if a{z, z) = 1 then U{z, z) = zz. We 
adopt the following general form of U: 

U{Z,Z) = - J2 TmnZ'^Z^ (15) 
m,n>l 



The expansion of $ around z = can be obtained from (0) in a direct way. The 
result is: 

$+(2, z) = -U{z, z)+vo + 27^e ^ t^/' (16) 

A:>0 

Here the non-harmonic term, U, is determined by the background charge, as is explained 
above, and does not depend on the form of the domain. It is the harmonic contribution 
that takes the boundary into account. The coefficients are 

tfc = ^j^^-'5,?7rfz {k>0) (17) 

They provide a complimentary (to Vk) set of parameters characterizing the domain D^. 
Here and below we use the same notation for the moments (to, tk, etc) as in Sec. 1. In 
the case cr = 1 they are equal to harmonic moments of the exterior domain D_ . Finally, 



^^0 



2 f 

— / log \z\ (r(z, z) d z. 

71 JD+ 



Note that the moments Vk (|T3D and to (0) also admit a similar contour integral repre- 
sentation: 

Vk = { z'^d^U dz {k>0), to = ^(fd,Udz (19) 

Z711 J-y Z7n J7 

which is a simple consequence of the Green formula. 

Since the potential and the electric field are continuous at the boundary of the domain, 
the two complementary sets of moments are related by the conditions 

= = 9^$" , z e-f (20) 

The generalized inverse potential problem is to find the curve 7 given the function a 
and one of the functions or $~, i.e., one of the infinite sets of moments. We will 
choose as independent variables the total charge vrto and the coefficients tk, tk [k > 1) of 
the multipole expansion of $ around the origin. Under certain conditions, which we do 
not discuss here, they uniquely determine the form of the curve, at least if the deformed 
curve to be described by tk is close enough to a fixed given curve. (Some important 
results concerning uniqueness of solution to inverse potential problems can be found in 



r^.) One such condition, specific for the case of non-uniform density of the background 
charge, is worth mentioning. For the problem to be well-posed, one should require that 
a{z, z) 7^ at any point of the curve 7. Otherwise a small perturbation of the curve 
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around the zero of a does not result, at least in the first order, in the variation of $. To 
avoid this difficulty, one may assume, for example, that a{z^ z) > 0. 

In other words, under the conditions mentioned above, {tk}'kLo is a good set of local 
coordinates in the space of analytic curves. The moments vj are also determined by 
{tfcjfcLo' moment of the interior can be regarded as a function of tk- For simplicity, 

we assume in this paper that only a finite number of tk and T^n are non-zero. In this 
case the series (|16D is a polynomial in z, z and, therefore, it defines the function for 
all z G -D+. As soon as a is real, to, are real numbers while all other moments are in 
general complex quantities. 

We now pass to the electrostatic variational principle mentioned in the beginning of 
this section. Consider the energy functional describing a charge with a density pi^z^ z) 
in the background potential $ (pA] ) generated by the charge distributed in with the 
density a: 

E{p) = -^ J Jd'^z dh' p{z, z) log |z - z!\p{z!, z) J d^z p{z, z) $(2, z). (21) 

The first term is the 2D "Coulomb" energy of the charge while the second one is the 
energy of interaction with the background charge. Clearly, the distribution of the charge 
neutralizing the background charge gives minimum to the functional. We denote density 
of this distribution by po^ Po = inside the domain and po = outside, so that the total 
charge density be zero. At the minimum the functional is equal to minus electrostatic 
energy E of the background charge minp£^(p) = — i?, where 

E = -\ 1 1 d^zd^z' a{z, z) log \z - z'\a{z', z) (22) 

'R J JD+ 

Varying p in (plf ) and then setting p = pq, we obtain eq. ([T6|). 

Let us find derivatives of E with respect to tj.. This can be done in different ways. 
One of them would be to take the derivative of ( p2D directly. In doing so, one should 
somehow control the contribution coming from an infinitesimal change of the domain 
under the variation of tk- To avoid this problem, we use the same argument as in the 
case (7 = 1 (see 0). Namely, instead of taking the derivative of (P^), we differentiate 
E represented as —£{p) in the form ( pi]) at the extremum. There are contributions of 
two kinds: one from the explicit t^-dependence of and another from the implicit 
tfc-dependence of pq: 



dE d£{p) 



dtk dtk 



-I d'za{z,z)dt,^Hz,z)- [d'z^^^ 
71 Jd+ J 6p{z, z) 



P=PO 



Spoiz,z) 
Stk 

P=PO 



Since we are at the extremum, the variational derivative in the second term is zero, and 
this term does not contribute. In other words, the derivative is equal to the partial deriva- 
tive of —S calculated at the extremum. Let us treat, for a while, Vq and tk as indepen- 
dent variables, then the partial derivatives of are especially simple: dt^^~^{z, z) = z^, 
9^(,$"'"(z, z) = 1. Plugging them into the integral, we arrive at the relations 

dE dE _ dE 

^ =^k, 7^ = Vk, = -to, (23) 

dtk dtk ovq 
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where the partial derivative with respect to is taken at fixed vq and tj {j ^ 0,k). 
Therefore, the differential dE is dE = J2k>oivkdtk + Vkdik) — todvQ. Note that there 
are two contributions to the t^-derivative of (p2D: one from (9(^$ and another from the 
change of the domain. The former is easily calculated to be ^Vk, thus the latter is ^Vk, 
too. However, it is not so easy to obtain this result directly (see the appendix for an 
example of such a calculation). A direct proof of the formulas for first derivatives in the 
case (7 = 1 can be found in 0]. 

It is more natural, however, to treat the total charge rather than fo as an independent 
variable, i.e., to apply the variational principle at a fixed total charge. This is achieved 
via the Legendre transformation. Let us introduce the function F = E + toVo, whose 
differential is dF = J^kyoi'^kdtk + Vkdik) + foc^^o- The integral representation of this 
function follows from (^): 

F = -\ n d^zd^z' a{z, z) log \z-^ - z'~^\a{z\ z') (24) 



TT^ J JD+ 

The function F plays a major role in what follows. It is a real-valued function of the 
moments to, ti, ^2, • • •• Rewriting (p3D in the new variables, we get the main property of 
the function F: 

dF dF _ dF 

otk dtk dto 

where the derivative with respect to tk is taken at fixed tj (j 7^ k). The very existence 
of the common potential function for moments implies the symmetry relations for their 
derivatives, dt„Vk = dt^Vn, dt„Vk = dt^Vn, which were first obtained in |I| (for the case 
0" = 1) within a different approach. 

To represent the result (^) in a more compact form, let us modify the potential $ 
and introduce $(2;, z) = ^{z, z) + 2to log \z\ + fo, so that 

z) = -- f d\ a(C, log - C"' (26) 



TT JD+ 

and 

F = — I d^z a{z, z) l'+(2, z) (27) 
27r JD+ 

The modification amounts to adding the neutralizing point-like charge at the origin. The 
expansion of the modified potential around the origin is 

<l+(^, z) = -Uiz, z) + 2to log 1^1 + itkz' + hz") (28) 

k>l 

From now on, this formula will serve as a definition of the variables tk- Likewise the 
function $ and its first derivatives are continuous on the boundary 7 = dD^. For z in 
the function $(2;, z) is harmonic. In terms of $, formulas (^) can be compactly 
written in the form 

V{z,z)F = ^-{z,z) (29) 

where we use the operator V introduced in (^), and z is assumed to be in D_. Both sides 
are to be understood as series in z and z. It is our assumption that they converge in D_. 
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To conclude this section, let us discuss homogeneity properties of the function F. 
Under the assumption that only a finite number of t^s are non-zero, we can substitute 
(P^) into and perform the term-wise integration. We get the following relation: 



2F = --^ / d^zUAU + tovo + y^itkVk + tkVk) , (30) 



where A = 4:dzdz is the Laplace operator. Generally speaking, (|30D does not have any 
definite homogeneity properties as a function of tk- However, for homogeneous functions 
cr, F is quasihomogeneous. For instance, if a{z,z) = with an integer M, then 

(PPP can be brought into the form 

4MF = -tl + 2Mtodt,F + ^(2M - k)itkdt,F + hdi.F) 

k>l 



where the relations ( P^D are taken into account. This relation reflects homoheneity prop- 
erties of the moments under the rescaling z ^ Xz and means that F is a quasihomoge- 
neous function of degree 4M. 

A more fundamental homogeneity property, valid for any cr, can be revealed by ex- 
tending the definition of the function F. Let us treat parameters T^^ defining U or a 
(see (|T5D) on equal footing with the variables t^. To wit, let us identify = T^o, ik = Tq^, 
to = Too, and allow F to be a function of all Tm„, m,n >0, regarded as an extended set 
of independent variables. Using the same variational argument as before, one can deduce 
the relations 



dF 1 



z''z"'a{z,z)d'z 



dTmn TT Jd+ 

which generalize (|25D. Then ( pOD acquires the form of the pure homogeneity condition: 

dF 

2F= Y: Tmnj^ (31) 

n,m>0 



The r-function can be defined as r = e^. This notation, however, is introduced here 
only for the purpose to stress the relation to Hirota equations. It is F, i.e., the logarithm 
of the r-function, which we only deal with in the sequel. Recall that the r-function 
itself does not have a good dispersionless limit. Only its logarithm, multiplied by a small 
dispersion parameter, makes sense in this limit. The bilinear Hirota equations for r are 
rewritten as highly non-linear equations for logr. In the next section, we derive these 
equations starting from the definition of the function F . 



3 Dirichlet problem and dispersionless Hirota equa- 
tions for the function F 

The dispersionless Hirota equations are certain relations between second derivatives of 
logarithm of the r-function with respect to tk- The second derivatives of F are coefficients 
of the expressions like D{zi)D{z2)F , where the operator D{z) is introduced in (^. Recall 
that we already know first derivatives of F, see (^9]). Therefore, what we need to find is 
V{zi, Zi)^{z2, Z2) for zi, Z2 G -D_. For this we use the following general argument. 
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Consider a small deformation of the domain obtained from the original one by 
adding to it an infinitesimal bump (of arbitrary form) with area e located at a point 
^ G 7. By denote variation of any quantity under such deformation in the first 

order in e. The potential generated by the deformed domain is then $ + Se{^)^, where 

5,(5) <l(^, z) = - log \z-' - rT (32) 

is the potential generated by the bump, as it is clear from the explicit formula (|26|) . 
Comparing the expansion of 6^(^^-^^{z, z) in z as z —>■ with (0), we find variations of 
under the deformation: 

Given any function A of the moments, its variation, Se{^)A, is 

k>l 

which can be written as 

5,(5)A = -a(e,eXe,0^, (33) 

So we see that for on the boundary curve 7 = ODj^ the operator 'D{^,S,) has a clear 
geometrical meaning: the result of its action on any quantity is proportional to the 
variation of this quantity under attaching the bump at the point ^. To put it differently, 
we can say that the boundary value of the function V{z, z)A is equal to 'k5^[z)A/ {ea{z, ^)), 
2; G 7. For functions A such that the series D{z)A converges everywhere in D_ up to 
the boundary, this remark gives a usable method to find the function ©(-z, z)A. To wit, 
this function is harmonic in D_, and its boundary value is given by (pSf). It is then just 
the subject of the Dirichlet boundary problem to find the harmonic function given its 
boundary value. This method appears to be especially useful when one is able to find 
the left hand side of (|33[) independently. 

As a very simple example, one can easily find 5e{£,)F = ^(y{^,^)^{^,^) from the 
integral representation of In this case, the harmonic continuation does not require 
any additional care since the function $ is already harmonic in D_. Whence we reproduce 
the result 



Let us turn to second derivatives. As it follows from (p3D , the boundary value of 
the function V{^,C,)'D{z,z)F = V{C,,^)^{z,z), regarded as a function of ^, is given by 
eq. (p2D for z in D_: 

V{C, l)V{z, z)F = -2 log \z-'-r'\, ^ e 7 (34) 

Here we do not discuss convergence of the series D{^)^ in D^, which is the necessary 
assumption for the harmonic continuation rule explained above. To be sure that the class 
of domains having this property is not empty, one may think of domains and background 
charges that are close enough to the uniformly charged disk, in which case the proof of the 
convergence is the matter of some routine estimates. The right hand side of (0) is not 
harmonic in D_ as it stands because of the logarithmic singularity at ^ = 2;. To find the 
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harmonic continuation of this function, we add to ( |3^ a function harmonic everywhere 
in D_ but at the point z, where it has the logarithmic singularity +2 log \ z^^ —^^^\ (thus 
cancelling the singularity of (|5^)), and equal to zero on the boundary. By definition of the 
Green function of the Dirichlet problem (see Sec. 1), it is the Green function G{z,^) (@). 
Hence we obtain formula (|^). Tending Z2 to oo and then separating the holomorphic part 
in Zl, one arrives at the expression (|^) for the conformal map. Separating holomorphic 
and antiholomorphic parts of (0) in both variables, and using (^, we obtain: 



log "''^'''"'^-' ^-kf + g(^OCfe)f (35) 

Zl — Z2 / 



log ( 1 - , I , A = -D{zi)D{z2)F (36) 

V w[Zi)w[Z2)j 

The limiting case of ( p5D as 2:2 ^ oo gives another useful formula for the conformal map: 

nj{z) = e-^^'o^ {z - {dt, + D{z))dt,F) (37) 

All these formulas have the same form for any function a. 

It is now straightforward to see how the Hirota equations come into play. They are 
obtained by plugging (||) into equalities of the type (^), (|3^) , thus excluding w{z) from 
them. In this way one obtains eqs. (§), (P), in which one recognizes dispersionless limit 
of bilinear Hirota equations for the 2D Toda chain hierarchy. Another Hirota equation 
does not even require (P) for its derivation: take three copies of eq. (|35|) for each pair of 
the points zi,Z2,z^, take the exp-function of both sides of each equation and then sum 
them up. One arrives at the equation 

{Zi - ^2)e^(^^)^(^^)^ + {Z2 - ^3)e^(-2)^(-3)F ^ _ ^^yDi.,)D(z,)F ^ q (gg) 

This is the dispersionless KP hierarchy (which is a holomorphic, with respect to the 
variables t^, sector of the Toda chain) written in the most symmetric form. Another 
version, known in the literature [1^, |Tl|, is obtained from here as z^ 00. 

The three equations, (H), (||), and (|38D (along with their bar- versions) are to be 
understood as an infinite set of relations for second derivatives of F. Altogether, they 
form the integrable hierarchy. These relations are obtained by expanding both sides in 
a power series in Zi and comparing the coefficients. For example, comparing leading 
terms in both sides of a.s zi, Z2 ^ 00, one obtaines the first equation of the hierarchy. 
It is the dispersionless version of the Toda equation: d^^i^F = e^*o^. (In terms of the 
field = f = dtgF, it is d"^ ^ (j) = (9t(,e^'o<^.) Next-to-leading terms give the equation 
dlt.F = KuFdl,^F. 

An important comment is in order. All the formulas for second derivatives of F, 
obtained in this section, are of the same form for any function a, i.e., a does not enter 
them explicitly. Indeed, consider, for instance, (|35|) , (|36|) . Their left hand sides, being 



expressed through the conformal map only, are determined mearly by the shape of the 
domain. For a given domain, they are the same for any density of the background charge. 
Clearly, this also holds true for the generating formula for second derivatives (^. Note 
that F itself does depend on the choice of a. 
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This remark can be reformulated as a covariance property of the second derivatives. 
Let F, ik be the dependent and independent variables related to another density function, 
a{z, z). As it follows from the above comment, second derivatives of F are covariant in 
the following sense: 

^ = ^ (39) 
dijdik dtjdtk ^ ' 

when the both sides are calculated for the same domain D^. A similar formula for 
mixed tj, tfc-derivatives also holds true. We stress that this covariance takes place only 
for second derivatives. For instance, as it is seen from (^), in general d^.F ^ dtjF, since 

the potential does depend on a for a given domain. 

The Hirota equations are invariant with respect to varying the density a. At the same 
time, the particular solution describing the conformal maps and the Dirichlet problem 
is determined by the background charge. We conclude that the function a parametrizes 
different solutions of the Hirota equations. We come back to this point again, within a 
different approach, in Sec. 5. Regarded as a function of all Tmn (see the end of Sec. 2), F 



solves the extended Toda chain hierarchy introduced in [13] in the context of the 2-matrix 
model. 

A short comment on the bilinear Hirota equations and their dispersionless limit is 
in order. Let ^ be a parameter. The equations are written for a r-function r/j which 
depends on the times and on the parameter h. The list of equations is parallel to the 
one for the dispersionless case. We use the same notation. Eqs. (||), (|) are replaced by 

= {z, - z,) (e-'^(^(-)+^(-))r,) {e^'^^m) (40) 

= (ziZ2)-'(e-'^(^*o+^{^i))r,)(e^(^*o+5(22))^^) (41) 

respectively. The dispersionful analog of eq. (|38D is 

[zi - Z2) (e-^(^(^i)+^(^2))Tn) (e-^^(^^Vn) + [cychc per-s of zi, z^, z^\ = (42) 

The operators e^^^^^ are shift operators, so the above equations are difference ones. To 
make this explicit, introduce the vector fields dx^ = —D{zi). In the variables Xi the 
equations aquire a more familiar form. For example, eq. (^) reads 

{zi - Z2)Tn{xi + h, X2 + h, X3)Tfi{xi, X2, X3 + k) + [cyclic per-s of 1, 2, 3] = (43) 

These bilinear difference equations for the r-function enjoy many remarkable properties. 
Eq. (^2]) first appeared in Hirota's paper [jl5|. For a review of the difference Hirota 



equations see e.g. |]TB| . 

As is clear from these formulas, the parameter h plays the role of lattice spacing. 
In the dispersionless limit, difference equations become differential ones. The limit is 
well-defined on the class of solutions such that there exists a finite limit 

F = lim^^logra (44) 
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Rewriting the Hirota equations for logr/i and extracting the leading terms ash —> 0, one 
arrives at the dispersionless Hirota equations for F. Different aspects of the mathematical 
theory of dispersionless integrable equations were discussed in |T^, |TB|, |TT], |T^. 



4 Integrable structure of conformal maps in the Ta- 
kasaki-Takebe formalism 

The approach of the previuous section is maybe the easiest way to see the link between 
Dirichlet problem and Hirota equations. However, such more customary attributes of 
integrability like Lax representation remain obscure, from this point of view. In this 
section we discuss an alternative approach which makes the integrable structure explicit 
via the Lax-type equations for the inverse conformal map, z{w). In our normalization, 
the general form of the inverse map is 

oo 

z{w) = + ^ UjW~^ . (45) 

j=0 

where r is real positive. For w on the unit circle, eq.(^) gives a parametrization of the 
curve. 

We introduce one more new object. Let S{z) be the analytic continuation of the 
function dzU away from the curve: 

S{z)=dM{z,z), ze'j (46) 

For analytic curves, S{z) can be proved to be holomorphic in some strip-like neighbour- 
hood of the curve 7. This function turns out to be a very useful object from technical 
point of view. We call it generalized Schwarz function, or simply Schwarz function for 
brevity. (The latter name is commonly used in the case a = 1, see e.g. the book ||19|].) 



Comparing its definition with the definition of the moments through contour integrals 
(0)) (0)' o'^^ obtains the Laurent series representation of the Schwarz function: 

00 J. 00 

Siz) = ^ ktkz'-' + - + ^kZ-"-'. (47) 

k=l ^ k=l 

We thus see that S{z) is the generating function of all moments tk, Vk- We regard it as 
a function of the moments tk'. S{z) = S{tk,z). Given the coefficients Tmn ( |15D and the 
Schwarz function, the curve 7 is determined by eq. (|i6|). 

Following the lines of ref. [Tl|], we introduce the dispersionless 2D Toda hierarchy 



starting with the following data. Suppose we are given with four Laurent series of the 
form 

00 00 00 

L = rw + J2 > M = Y, ktkL^ + + I] VkL'^ (48) 

j=0 k=l k=l 



00 00 



j=0 k=l k=l 
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The coefficients tQ,tk, tk are regarded as independent variables while all other coefficients 
(r, Uj, Uj, Vj, Vj) are dependent variables. Imposing some conditions on the Laurent series, 
we are going to study the latter as functions of the former. 

Recall now the Takasaki-Takebe theorem |]ll|. Let the Poisson bracket {, } be defined 

as 

and let f,ghe a. canonical pair, {f{w, to), g{u), to)} = f{w, to). Suppose the four Laurent 
series L, L, M, M introduced above are connected by the functional relations 

L = r\L,M), M = g{L,M) (51) 

Then the following is true: 

(A) The pairs L, M and (Z)~\ M are canonical: {L, M} = L, {L-^, M} = L-^; 

(B) The four Laurent series satisfy the Lax-Sato equations 

where X stands for any one of L, L, M, M and 

if„= (L")>i + ^(L")o, ^„ = (L")<_i + l(L")o, n>l (53) 

(C) There exists a function F of t^, ik, and to such that Vk = dt^F, Vk = dt^F, and the 
function F obeys the dispersionless Hirota equations. 



This theorem is proved in [|Tl|]. The proof of (A) consists in differentiating (|5lD with 



respect to if, to and combining the results. Taking derivatives of (^) with respect to 
tn,in and using (A), we obtain (B), where (...)>i ((...)<_i) means the part of the Laurent 
series with strictly positive (strictly negative) powers of w, and the constant term is 
denoted by (...)o. To prove (C), one should take dM/dtk at fixed w. Using the Lax-Sato 
equations (|52D , one gets dt-V^ = {wL^dHj)^. Further, this representation, together with 
(p3|), implies the symmetry dt v^ = dt^Vj, which, in its turn, implies the existence of a 
common potential function for Vk- The proof of the Hirota equations for F from these 



data is implicit in |TT|. More explicit arguments are given in (see also |0|). 

It then follows that L and L are the two Lax functions of the dispersionless Toda 
hierarchy. Equations of the hierarchy are written for the "potentials" r, {«„} 
as functions of the times to, {t„}, {t„}. These times appear expicitly in the Laurent 
series M, M which are the Orlov-Shulman functions ||20| of the hierarchy. The Takasaki- 
Takebe theorem says that solutions to the hierarchy are in one-to-one correspondence 
with canonical pairs of functions /, g. 

Remarkably, this formalism turned out to be very convenient for describing the inte- 
grable structure of conformal maps. The identification of the objects related to conformal 
maps with the above Laurent series goes as follows: 

L{w) = z{w) , L{w) = z{w''^) , M(L) = zS{z) , M(L) = zS{z) 
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(for a series f{z) = J2 fj^^ we write f{z) = J2 fj^'')- Together with S{z) = d^U^z, z) this 
means that we have the following functional relations between L, L, M, M: 

M^L?^, M^L?^ (54) 

Notice that (|54D gives an equivalent description of the canonical pair from the Ta- 
kasaki-Takebe theorem. Namely, (5i) tells us that U{L,L) is generating function of 



the canonical transformation L,M — > L ,M and, therefore, implicitly determines the 
corresponding canonical /, (/-pair. Indeed, let in general M = F{L,L), M = G{L,L), 
and find the Poisson bracket between and M. For brevity we write L = ip{L,M), 
where is a function implicitly determined by the equality M = F{L, L). We have: 

dLdM dL dM _ dipdG dip dcp f dG dG dip\ _ dcp dG 
dZdM ~ dM~dL ~ dZ'dLdM ~ dM \dL ^ 'dLdZJ ~ "dU'dL 

Therefore, 

where the last equality follows from dip/dM = [dL / dM) i^const = {OF / dL)~^ . Now, on 
taking F,G from (0), we get {L~^i M} = L~^i i.e. the transformation (Q) is indeed 
canonical. 

The Poisson bracket {L, dLU{L, L)} = 1 follows from {L, M} = L. This can be 
rewritten as 

{L, L} = — ^ (55) 
o-(L, L) 

which is one of the forms of the so-called string equation. 



5 Conclusion 



We have demonstrated that some classical problems of complex analysis have an inte- 
grable structure which is the same as the one arising in the dispersionless limit of soliton 
equations. Given a simply connected domain in the complex plane and an auxiliary 
real-analytic function a{z,z) defined all over the complex plane, we have constructed a 
solution to the dispersionless 2D Toda chain hierarchy. The independent flows of the 
hierarchy are moments of the domain, or multipole coefficients, defined with respect to 
the function a which in the context of potential theory plays the role of density of back- 
ground charge. The set of dependent variables is then the complimentary set of multipole 
coefficients. The energy of the background charge in the domain, regarded as a function 
of the moments, is (logarithm of) the r-function of the hierarchy. This function obeys the 
dispersionless Hirota equations. Taking derivatives of this function with respect to the 
moments, one obtains formal solutions to the conformal mapping problem, the Dirichlet 
boundary problem, and to the 2D inverse potential problem. An interesting question, to 
be discussed elsewhere, is an inverse problem connected with the scheme outlined above: 
given a solution to the dispersionless Hirota equations, associate with it a boundary or 
conformal mapping problem, and to restore the function a. Such an interpretation may 
be helpful for better understanding dispersionless integrable hierarchies. 
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Appendix 

In the appendix we give a direct proof of the formulas (^) for first derivatives of the 
function F with respect to t^. We use the notation of Sec. 2. 

Consider the variation 6kF {k > 0) of the function 

F=—[ ^{z,z)a{z,z)d'^z 
271 Jd+ 

under a small change —>■ tk + Stk'- F ^ F + S^F. All other moments tj with j ^ k, 
j >0, are kept constant. We have: 

6kF = ^f 6k^+{z,z)a{z,z)(fz + ^ f ^{z, z)a{z, z)d^z 

Z7l JD+ 2,71 J&kD+ 

where 5^-0+ is the variation of the domain due to the change of t^'s. (To be more precise, 
/(5feD+ stands for the difference Jn+itk+stk) ~ In+itk)-) -^i -^2 we denote the first and 
the second terms of this formula. Ji is easy to calculate. Since Sk^~^{z, z) = z^5tk + z^5tk, 
the first term is 

h = — 5k^^{z, z)cr{z, z)(fz = -{vkStk + VkStk) 
Ztc JD+ 2 

Let us turn to the second one, h = ^ hkD+ ^{^^ ^)a{z, z)(P'z. Since the integral goes 
over a small neighbourhood of the curve 7 = i9-D+, we can substitute $ by $~, which, by 
virtue of (pUD, coincides in this neighbourhood with up to second order terms. Recall 
that ^~{z, z) = Vq + 27leuj{z), where the function 

k>l ^ 

is holomorphic in and regular at infinity. Under our assumptions, all singularities of 
uj are in a domain B C Z)+ such that R = -D+\-B is a strip-like domain. In other words, 
uj{z) can be analytically extended across 7 to the strip-like domain R. Moreover, recall 
that we are working in the space of analytic curves. It can be proved that for curves that 
belong to a small but finite neighbourhood of the point in this space that corresponds to 
a given curve, the domain B can be chosen to be independent of t^. 

We can write, in the first order in 61^: 

l2 = — [ a-{z,z)(fz-\ — Tie [ uj{z)a{z, z)d'^z 

27[ J5kD+ 71 JSkR 
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The first integral in the right hand side is the variation of Iq, and, therefore, equals zero 
as soon as to is kept constant. So, we are left with the second one, 

I2 = —Tie [ uj(z)a(z,z)(fz Tie [ Lj(z)a(z,z)(fz = 

TT J D+{%+6t,,)\B TT JD+{t^)\B 

r f dz 

= Tie i uj{z)dJJ{z,z)—Tlef uj{z)dJJ{z,z) 

JdR{tk ) ZTTZ JdR{tk)2m 

The second line follows from the Stokes formula. The boundary of R consists of two 
curves: inner and outer ones. The inner boundary is the same in both integrals, so its 
contribution cancel. The outer boundary is 7(tfc + 6tk) in the first integral and 'y{tk) in 
the second one. Therefore, using the definition of the generalized Schwarz function (P6|), 
we can write: 

r dz f dz 

l2 = Tlei -—uj{z)S(tk + Stk,z) -Tie f -— Lj{z)S(tk, z) 

Since S{z) is holomorphic in R, the contour can be taken to be 7 = '~f{tk) in both integrals. 
Thus we obtain, by virtue of (^Tj): 

^ /r f dz . .dS(z)\ 1, ^ _ \ 
I2 = Tie \dtk <J> — uj{z)— — I = -{vkdtk + Vkdtk) 

Summing up the contributions Ji and I2, we get 6kF = VkStk + VkSik, which is the desired 
result for k > 1. The calculation of the to-derivative can be done in a similar way. 
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